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REFINED EIGENVALUE BOUNDS ON THE DIRICHLET 
FRACTIONAL LAPLACIAN 

SELMA YILDIRIM YOLCU AND TURKAY YOLCU 


Abstract. The purpose of this article is to establish new lower bounds for 
the sums of powers of eigenvalues of the Dirichlet fractional Laplacian operator 
(—A) Q / 2 |f 2 restricted to a bounded domain SI C K d with d = 2, 1 < a < 2 
and d > 3, 0 < a < 2. Our main result yields a sharper lower bound, in the 
sense of Weyl asymptotics, for the Berezin-Li-Yau type inequality improving 
the previous result in [36]. Furthermore, we give a result improving the bounds 
for analogous elliptic operators in [19]. 


1. INTRODUCTION 


Fractional Laplacian operators are usually considered as the prototype of the 
non-local operators [11] and have recently garnered much attention in many appli¬ 
cations in mathematics and physics. Related problems usually lie at the interface of 
probability, stochastic processes, partial differential equations and spectral theory 
such as [5, 11, 12]. For some applications, we refer the reader to graphene models 
[15], obstacle problems [29], non-local minimal surfaces [9]. 

In this article, we study eigenvalues of the fractional Laplacian operator (—A) a / 2 
defined by 

(-A) Q/2 fa = Aj a) fa in fl, 
fa =0 in R d \fl 


where Q is a bounded connected domain with smooth boundary in R d , for either 
d = 2, l<a<2ord>3 and 0 < a < 2. Since is bounded, the spectrum of the 
fractional Laplacian is discrete and eigenvalues (including multiplicities) 

can be sorted in an increasing order. 

Unlike Laplacian, fractional Laplacian is a non-local operator such that for suit¬ 
able test functions, including all functions in it € Gg°(R d ), it is defined as 


(-A ) a/2 u(x) = A d , a lim 
£-> 0 + 


i(x + y) — u(x) 


'{\y\>c} 


I y\ 


d-\-<y. 


dy , 


( 2 ) 


where Ad,a is a well-known positive normalizing constant. From a probabilistic 
point of view, (—A)“/ 2 can be considered as the inhnitesimal generator of the semi¬ 
group of the symmetric a—stable process, denoted by A' t , with the characteristic 
function 


e -t|Ml“ = E ( e ivx*) = f e ^ y p[ a) {y)dy, t > 0, (3) 

JR d 
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where p[ a \x, y) = p[ a \x — y) is called the transition density of the stable process 
(or the heat kernel of the fractional Laplacian). While explicit formulae for the 
transition density of symmetric a-stable processes are only available for the Cauchy 
process (a = 1) and the Brownian motion (a = 2), these processes share many of 
the basic properties of the Brownian motion. Another process of importance is 
the Holtsmark distribution (a = 3/2) which is used to model gravitational fields 
of stars (See e.g., [38]). Stable processes do not have continuous paths which is 
related to non-locality of the fractional Laplacian operator [5, 8]. 

When fractional Laplacians involved, some of the known methods fail because 
of the fractional power and non-locality of such operators. This drawback can be 
evaded by using the Fourier transform definition. Recall that the Fourier transform 
and its inverse are defined as 

F[u]{p) = u{n) =c d [ e~' lfJ, ' z u(z) dz, ,F - 1 [u](z) = c d f e lfl ' z u(p) dp, 
jR d J R d 


where c d = (27t) = is the normalizing constant. Interestingly, the fractional Lapla¬ 
cian operator on fl C can be defined as a pseudo-differential operator as 

(-A)“/ 2 | nW = F~ x M a f[u}] , 0 < a < 2, ueH% /2 (n). (4) 

Here, H^’ 2 (Q) denotes the Sobolev space of order a/2. When = R d , one can 
look at Proposition 3.3. [13] for the proof of the equivalence between the definitions 
in (2) and (4). 

There is an extensive literature devoted to the inequalities involving the eigen¬ 
values of the Dirichlet Laplacian operator, which can be regarded as the fractional 
Laplacian when a = 2. One may consult the articles [2, 3, 4, 14, 17, 20, 24, 25, 30] 
and references therein for a through literature review. It is worth pausing here for 
a moment to consider the Dirichlet Laplacian results relevant to our main result. 
The first such result is the Li-Yau inequality that provides a lower bound for the 
sums of eigenvalues, sharp in the sense of Weyl asymptotics [26]. The authors gen¬ 
eralized this result in [33] by obtaining the following Li-Yau type inequality for the 
eigenvalues of the Diriclet fractional Laplacian operator: 


k 


3 =1 


(a) 


> (47r) 



- 1+7 


( 5 ) 


where |fl| represents the volume of f l and r(:c) denotes the Gamma function T(a;) = 
J^°t x ~ 1 e~ t dt for x > 0. One may also look at [16] for a Li-Yau type inequality 
involving the eigenvalues of the (massless) Klein-Gordon square root operators 
(—A) 1 / 2 |n, (i.e., the case a = 1). To look at this inequality from a different 
perspective, one may take the Legendre transform of the following result by Laptev 
[23] and obtain (5): 


^( 2 -A;.“ ) ) + <(4^-i 


a M ..1+4 

a + dr(l + f) 


( 6 ) 


When we set a = 2 in (6), we recover an earlier result by Berezin [6], which supplies 
the Li-Yau inequality after an application of the Legendre transform. Thus, in what 
follows, we call (5) as the Berezin-Li-Yau inequality. 
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In [27], Melas proved the following improvement to the Berezin-Li-Yau inequality 
(a = 2): 


£A?> >4tt 

3 =1 


r (l+ I) 


2 + d 


l«l 


fc 1+ 3 + 


| 0 | 


24(2 + d)l(n) ’ 


( 7 ) 


where 1(12), the moment of inertia, is defined by 


I(f]) = min / \z — y\ 2 dz. 
ve«. d Jn 

By a translation of the origin and a rotation of axes if necessary, in the sequel, we 
assume that the origin is the center of mass of Q and that 


Z(ft)= f \z\ 2 dz. (8) 

Jn 

Melas type bounds and their many variants and extensions have recently attracted 
a lot of attention, see for instance [18, 21, 31, 32, 33, 34, 35, 36, 37]. In particular, 
in [33], the authors obtained a refinement of (5), stating that 


£ A$ a) > (4 t r)t -A- ( ' k 1+ 

fri J ~ a + d \ |n| J 

a |nr-^r(i + f)^- ,, 1+ ^ 


48 (a + d) (4 7 r) 1 -f 1(12) 




( 9 ) 


Remark that a = 2 in (9) recovers Melas’ bound in [27] for the Dirichlet Laplacian 
eigenvalues. Unfortunately, it is not easy to take Legendre transform of (9) to find 
a similar improved Berezin type bound in the case of the fractional Laplacian. 

For 0 < £ < 1, authors also proved in [36] that 


E (E 1 ) 


e 


> (47t)" 


a£ + d 


E0±|n d fci+ 


i«i 


16 (a£ + d) (47r) 1- ^ 1(12) 

[H| 2 -^r(i + |) : # u+ , 

„ ,v ft 


:l-2 

d 


6A0(a£ + d) (47t) 2 2 X(I2) 2 


( 10 ) 


The main purpose of this article is to obtain analogous but sharper bounds with 
Dirichlet Laplacian replaced by Dirichlet fractional Laplacian. The first step, in¬ 
spired by results in [31, 32, 33, 36, 37], is establishing the following result: 


Proposition 1. For k > 1, and either 1 < a < 2 and d=2 or 0<a<2 and 
d > 3, the eigenvalues {A) a ^}^E 1 of the fractional Laplacian operator (1) defined on 
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0 C satisfy 




0 =1 


g-1 

d\ d 


2(a + d) 

(47r)2-§ X(fi)3 

5a 

|ii| 1 -^r(i + 

16(a + d) 

(4 7 r) 1 -f I(fi) 

a 


16(a + d) 

(47r)3 _ t 1(0)1 


*-2 

d 


k 1+ —. 


(11) 


Note that the constants in the leading terms on the right side of (11), which is a 
fractional version of Weyl’s law, are optimal due to a classical result of Blumenthal 
and Getoor [7]. While the constant in the second term in (11) is still sub-optimal, 
the estimate in (11) is substantially stronger than previous known results in [33, 36, 
37]. In addition, we recover the lower bounds in the case of the Dirichlet Laplacian 
[32] when we set a = 2 in (11). 

In [12], Chen and Song obtained that 

y (12) 

for each j and any constant t £ (0,1]. Thus, Proposition 1 along with an application 
of (12) leads to our principal result: 


Theorem 1. For k > 1, 0 < £ < 1 and either 1 < a < 2 and d = 2 or 0 < a < 2 
and d > 3, the eigenvalues of the fractional Laplacian operator (1) defined 

on O C satisfy 


t(Af)'>(4^-yCf r( h +i) )'^ 


i =i 


at + d\ |0| 

at |S2|5-®Vt(1+ L) d kl+ «e=i 


2 (at + d) 

(47t)2 °2 1(0)2 

5 at 

|0| 1 -^P(1 + |)^ 

16 [at + d) 

(47r) 1- ^ 1(0) 

at 

o o , a£-3 

|0||- 2 ir 1 r(l + f) d 

16 (at + d) 

(47r)i-^T(0)l 


k 1+ — 


(13) 




On a side note, the proof of Proposition 1 consists of a very delicate application of 
Steffensen’s type inequalities, which is mainly about the comparison of the integrals 
on the subsets of interval [0, oo). 

The article is structured as follows: In Section 2, we revisit the relevant facts 
about the eigenvalues and eigenfunctions of the fractional Laplacian operator. After 
providing the proof of an auxilliary lemma that plays a crucial role in proving 
Proposition 1, we present the proof of our main results in Section 3. Finally, we 
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end Section 3 with a remark which extends the main result even further for some 
elliptic operators studied in [19, 28]. Please see Remark 2 for details. 

2. Preliminaries 

In this section, we give an overview of the definitions and tools that are essential 
to establish the estimates in (11). Even though, these were previously studied in 
[31, 33, 36, 37], we include them so that the article is self-contained. By using 
Plancherel’s theorem, one can show that the set of Fourier transforms {4>j}'jLi of 
{<l>j} < jLi forms an orthonormal set in L 2 (R d ) since the set of eigenfunctions 
is an orthonormal set in L 2 (S7). To ease the notation, we set 

k k ,, 2 

> o. 


®k{lA ■= E l4(X)| 2 = Ew E [ e lZ '' 1< t > o{ z ) dz 

i=i ^ W ’ l=i 

Because the support of tp 3 is f1, the integral is taken over 17 instead of ' 
changing the sum and integral and using 110^ 11 2 = 1, we derive 


/ "MlO dn = k. 

J R d 


(14) 
Inter- 

(15) 


Observe that 


E x j = E<^ EE) = E^i> (- A )“ /2 ^-) 


1=1 


3= 1 

k 


3 =1 


= ^(^,^- 1 [H“^i]]) = E [ H“|4'(m)I 2 ^ (16) 

i=i 3=1 jRd 

= [ \y\ a <5>k{v)dy,. 

J R d 


Application of Bessel’s inequality gives an upper bound for 


®k(v) < 


1 


(27r) a 


\e~iz-M\ dz = 


in 


Next, we find an estimate for |V$fc|. Note that 

k 




(27r) d 

:= n d . 

(17) 

| 2 dz - 

J(O) 

(27r) d 

(18) 


In view of Holder’s inequality and utilizing (17) and (18), we arrive at the following 
uniform bound: 


|V<Mm)|<2 (£|0j(m)I 
0=1 


1/2 / \ 1/2 

k 


E i v ^-(m)i 2 

i=i 

< 2 { 2 Tr)~ d y/\n\l(n) ■- p. 


(19) 


Let B r {z) := {y £ R d : \y — z\ < R} designate the ball of radius R centered at z in 
and oj d denotes the volume of d dimensional unit ball -Bi(z) in given by 


d 

7T 2 


U d = 


r(i + |)- 


( 20 ) 
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Now assume that R is such that |fi| = uidR d ■ That is, -Br(O) is the symmetric 
rearrangement of fl. Note that 


i(n) > 


■dz = 


>b r ( o) 


du) d 
d + 2 


R d+2 = 


d + 2 




roughly giving 

p> { 2 n)- d uf* . ( 21 ) 

Let denote the decreasing radial rearrangement of $fc(/z). There exists a real 

valued absolutely continuous function ipk ■ [0, oo) —)• [0, fid] such that 

= vM)- ( 22 ) 

By using Polya-Szego inequality, one can show that 

0 < < p. (23) 


For more details, see for example [33]. 


3. Proof of Proposition 1 


Before diving into the proof of the main results, we present the following surpris¬ 
ing sharper inequality which will be the main ingredient in the proof of the sharper 
lower bound in (11). Our method of proof has been previously explored in several 
articles, for instance [27, 31, 32, 34, 36], with crucial differences. 

Lemma 1. For either d = 2 and 1 < a < 2 or 3 < d £ N and 0 < a < 2, a > 0, 
b > 0, we have the following inequality 

a d+ a > ^±ff a d b a _ “ b d+a + a b d +a -3 £ 2a + ( a _ b f_ (24) 

d d d 

A direct but lengthy proof of this lemma is given in [36]. Here, we shall give a 
more intuitive and rigorous proof using convexity. 


Proof. First, let us show that 

h{x) := dx d+a — (d + a)x d + a — a(2x + l)(ai — l) 2 > 0. (25) 

for x > 0, either d = 2 and l<a<2or3<dsN and 0 < a < 2. 

Case 1: Assume that d > 3 and 0 < a < 2. Observe that h can be written as 
h(x) = x 2 g(x ) where 

g(x) = dx d+a ~ 2 — (d + a)x d ~ 2 — 2ax + 3a. 


Differentiating we get 

g\x) = d(d + a — 2)x d+a ~ 3 - 

g"{x) = x d ~ A d(d + a — 2){d + a — 3) 
Note that g"(xd, a ) = 0 where 


(d + a){d — 2)x d 3 — 2a, 

_ (d + a)(d — 2)(d — 3) \ 
v d(d T a — 2)(d -f- a — 3) J 


__ ( (d + a)(d — 2)(d — 3) \ 1/a 
d,a \d(d + a — 2)(d + a — 3) J 

When x > we have g"(x) > 0, implying that g is convex on [xd,a, oo). Thus, 

g(x) > g{ 1) + g\l){x - 1) = 0, 
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since g(l) = 0 and g'{ 1) = 0. That is, g(x) > 0 on \xd, a , oo). In particular, 
g(xd.a) > 0. On the other hand, when 0 < x < Xd, a , we have g"(x) < 0 yielding that 
g' is decreasing on [0,ird iQ ]. This implies that g'{x) < </(0) = —2a < 0, meaning 
that g is decreasing on [0,This leads to g(x) > g(xd, a ) > 0 for x € [0, Xd >a \- 
Hence, g{x) > 0 for x £ [0, oo). Therefore, we deduce that h{x) = x 2 g(x) > 0 for 
x £ [0, oo). 




Figure 1. Graphs of g{x) for d > 3 and d = 2, respectively. 


Case 2: Now, assume that d = 2 and 1 < a < 2. Then h becomes 

h{x) = 2cc 2+a — (2 + a)x 2 + a — a(2x + l)(x — l) 2 . 

As before, we can write h(x) = x 2 g{x) where g{x) = 2x a — 2ax + (2a — 2). Note 
that if a = 1, then g(x) = 0. Differenting again, we obtain 

</( x) = 2ax a ~ 1 — 2a and g”{x) = 2 a{a — l)a; a-2 . 

Notice that for 1 < a < 2, g"(x) > 0, meaning that g is convex. Since <?(1) = 0 
and g'( 1) = 0, g(x) > g( 1) + g'(l)(x — 1) = 0 implies that g(x) > 0 and, therefore, 
h(x) > 0. 

Setting x = a/b in (25), multiplying through by b d+a /d and rearranging the 
terms, we conclude (24), as desired. □ 

Remark 1. When 0 < a < 1 and d = 2, the inequality above fails to hold, therefore, 
we do not resolve this case in this manuscript. 

With Lemma 1 in hand, we are now ready to prove Proposition 1. 

Proof. Assume that (15)-(19) hold and d > 2. Consider the decreasing, absolutely 
continuous function <p k : [0, oo) —► [0, oo) defined by (22). We know that 0 < 
—ip' k (t) < P for t > 0 where /3 > 0 is given by (19). Since tp k { 0) > 0 due to (14) let 
us first define 

e ‘ (a) := (cir a ) ■ (26) 

Note that <d k is positive, 0^(0) = 1 and 0 < —Q k (a) < 1. To simplify the notation, 
we also set 9 k (a) := — 0' fc (a) for t > 0. Hence, 0 < 6 k (a) < 1 for t > 0 and 

/•OO 

/ 0 k {a) da = 0^(0) = 1. 

Jo 

Now, set 

/»00 POO 

Ck= a d_1 0fc(a)da and Vk — I a a+d_1 


o 


0 


0fc(a)da. (27) 
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Using (15) we get 

k= I dy = f dy. = doj d [ a d ~ 1 ipk(a)da. (28) 

J R d JR d Jo 

Moreover, since the map y i—> \y\ a is radial and increasing, by (16), we obtain that 
k 

= [ \v\ a 

3=1 jRd 

> [ imi 

JR d 

roo 

= dojd / a a+d ~ 1 (pk(a) da. 

Jo 

Substitution of (26) into (27) yields 

,d foo pd k 


Ck = 


Vk = 


Mo) d+1 Jo 

pa+d 

<Pk{ o) 


/*oo 

/ a d ~ 1 ipk{a) da = 

Jo 


du] d <Pk( 0 ) d+1: 


a+d +1 


, oa+d sr^k \(a) 

a° l+d - 1 ipk( a ) da < J_1 J 


dio d ipk{0 ) 


a-\-d -\-1 


Observe that application of Fubini’s theorem together with 

pOO 

&k(b) = / Skip) da 

Jb 


leads to 


_ t-\-d 


t -|- d , 


9k(a)da = J i^J b t+d 1 db'j 9k{a) da 

•OO \ 

9k{a) da ) db 


b t+d 


= / 6 *+<*- 1 © fc ( 6 ) rfb, 

Jo 

which together with y = 0 and y = a respectively yields 

pOO pOO 

/ a d 9k(a) da = (kd and / a a+d 9 k (a) da = ijk (ct + d). 
Jo Jo 


Notice that 


(a d - l) (9 k (a) - l[ 0 ,i](a)) >0, a £ [0,oo). 


(29) 


(30) 


(31) 

(32) 


Integrating (32) from 0 to oo gives 

r°° i 

J a d 9 k {a)da>-j-^= 7d(0), 

where 7 d ■ [ 0 , 00 ) —► ( 0 , 00 ) is defined by 

rx+l 

7 d(x) — / a d da. 

J X 

Since 7 ^ is continuous and non-decreasing and 7 ^( 2 ;) —> 00 as x —> 00 , the Interme¬ 
diate Value Theorem provides us with the existence of r > 0 such that 

pT+1 pOO 

a d 9k (a) da, 


/ r-t-i pc 

a d da = j 


a d da = 


0 
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which, by (31), concludes that 


r-T+l 


a d da = d(k- 


(33) 


Now consider the polynomial 

T(x) = x a+d - u lX d + is 2 = x d (x a - ui) + v 2 

where 


ot-\-d _ —OL-\-d 


Vl = 


(t + l) a+d - 
(t + V) d — T d 


(rr- _i_ 1 \a+<J _ ~-a+d 

>«. 


are chosen so that T(r) = 0 and T(t + 1) = 0 and T remains negative on (r, r + 1) 
and positive on [0, oo)\[t, t + 1]. It is immediate to observe that 


T{a) (l[ T , T+ i](a) - 9 k (a)) < 0 on [0,oo). 
Integration of (34) on [0, oo) leads to 


(34) 


/ r+l p oo / roc pr+1 \ 

a a+d da < J a a+d 9k(a) da — v\ yj a d 9 k (a)da — J a d daj, 


simplifying to 


/ r+l /»oo 

a a+d da < J a a+d 9 k (a) da. 


Using (31), we infer that 


/*T +1 

J a a+d da < r]k (a+ d). 


Observe that 


rr+l 


d(k = 


da> a d da = 


d+ 1 ’ 


(35) 


(36) 


(37) 


Notice that (24) gives the key inequality in the proof of this lemma. Indeed, inte¬ 
grating (24) in p from r to r + 1 we obtain 


/■r+l 


a «+ d da> ^±l h a 


rr+1 


a d da — — b a+d 
a 


/»T + 1 

+ — b a + d ~ 3 / (2a + 6)(a — b) 2 da. 

d Jr 


(38) 


Note that [32] 


FT -fTti 

/ (a — b) 2 da> min / (a — b) 2 da =—. 
Jr r>0, b>l/2_/ T 12 

/»r+l /»r+l 

/ a (a — b ) 2 da > min / a(a — b) 

Jr r>0,b>l/2Jr 


T + 1 i o 1 

— b) 2 da > -b 2 -6-1— 

“2 3 4 


and so, we have 


J (2a + b) (a — b) 2 da > b 2 — ^6 + ^ 


(39) 
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for any b > 1/2 and r > 0. Since (C kd ) 1/,d > (d + 1) x ! d >3 1//2 > 1/2 due to (37), 
setting b = (Cfcd) 1/,d and using (33) and (39), we deduce that (38) yields to 

1 


V k > 


j(C kd) 

5a 




d(a + d) 


(Ckd) 


i+- 


-(C kd) 


i+- 


:(Cfc d) 


(40) 


i+ £ 


4 d(cn + d) 2 d(a + d) 

Equations in (30) combined with (28) and (29) turn (40) into 

k 


E^> 

t=i 


a + d 
a 


u> d d <£fc(0) 2 k 1+d 




^fc(O) 1 


k l+l - 




(41) 


a 


P-3, 


r/3 d ^ fc (0) d <* 


2(a + d) 

To finish the proof of Proposition 1 we shall minimize the right side of (41) over 
(pk( 0). To prove (11), we show that the function, which we call d(x), on the right- 
hand side of (41) with x := <pk( 0) > 0 decreases on (0, fid]. By (17) we know that 
0 < x < f Id- First, define 

d( x) = $i(x) + d 2 (:r) 

where 

. , . dfc 1+ a _ a 

di(x) = -5- x d + 


(a + d) t 


y — 1 
d 


and 


d 2 (x) = 


a 


k l+L 


„3— - 


(a + d)/3u d c 
i 5 ak 1+ ~3~ 


‘2 — - 


(42) 


(43) 


2 (a + d)f3 3 uj d d + d)/3 2 cu d d 

Thus, it is enough to show that di and d 2 defined by (42) and (43) are also de¬ 
creasing on the interval (0, fid]. Differentiating di and d 2 > we observe that di(x) 
and d 2 {x) are decreasing on the interval (0,a,’i) and (0, * 2)1 respectively, where 


Xi = 


dpk d 


d 

d+1 


,2 (d — a + 1) 0 J d 


x 2 = 


5(2d + 2-a)/3fca 
, 2(3d + 3 — a)u) d 


d 

d+l 


Hence, we particularly obtain that d is decreasing on (0, fid] when we have fid < 
min {x -[, x 2 ) for any k > 1. Since x H > T(x) is increasing for x>2we obtain that 


r { 1+ 2 ) >r(2) = i. 


(44) 


In view of (44), (3 > (2n) d u> d d |fl| a and definition of uid given in (20) we observe 
that 


X\ > 


d|f2|^ fca 


> 


[r(l + f)p Ifll^fcE d+1 


, (d — a + 1) (2n) d ix d 


(27r) 


d+1 


>0*, 
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as k > 1 and d — a + 1 < 2d. Similarly, we obtain that 

(5(2d + 2 -a)\n\^k k A 41 ( [r(i + f)] f 

X2 ~ \ 2(3d + 3 — a) (27r) d wf / ~{ ^ d+1 


d 

d+1 


> fid, 


as k > 1 and 5(2d+2 — a) > (3d+3— a). In conclusion, we obtain that d(x) > d(£ld) 
as d is decreasing on (0, f2d]. Substitution of /3 = 2(27r) _d ^/|f2|X(fi) given in (19) 
together with (pk{ 0) = turns (41) into (11). □ 


Remark 2. In view of the recent work [19, 28], it is worth noting that one can easily 
extend this for elliptic operators £/ defined by a kernel / 


£ju{x) 

where / satisfies 


lim 


/ (u{x + y)-u{x))f(y)dy, 

(45) 

’{\v\>A 


f(y)>v^ d+a , 

(46) 


Ad,a is the normalizing constant in the fractional Laplacian definition (2) and er > 0. 
To this end, let us consider the eigenvalue problem defined by 


£J O'!y Ay l/)y 111 11, 

4>j =0 in R d \f2 


(47) 


It is shown in [28] that the spectrum of £f is also discrete and the eigenvalues 
(including multiplicities) can be sorted in an increasing order. Also, the 

set of Fourier transforms {bj}°8_:L of {bj}]pLi forms an orthonormal set in L 2 (R d ) 
since the set of eigenfunctions {bjlyki is an orthonormal set in L 2 (fl). Note that we 
use the same notation for eigenvalues and eigenfunctions to illuminate the striking 
similarities though they might be different for each £f. Defining T*, as in (14), we 
obtain (15) immediately. However, (16) needs to be re-written as the following 
inequality 

k k k k „ 

y = y "j ifiji^j^j) = y 'x &r = y i / Qai^) ^(aoi 

3 = 1 3 = 1 3=1 3 = 1 JWLd (48) 

> o- f |/r| a 4>/ c (/i) dfi. 

Jwi d 


where we used the fact (Proposition 3.3. in [13]) that 


Q»(ld)= (1-cos {yn))f(y)dy>aA d , t 

JwL d 


lK. d 


1 — cos(v ■ u) , . 

i g ,k: ^■ 


Having (48) in hand, (29) changes as follows: 

k » 

J2Xj>a \n\ a $ k (y,)dn 

3=1 jRd 

><j [ | n\ a <f>* k {n)diJ, 

JR d 

/»oo 

= adujd / a ot+d ~ 1 ipk{a) da. 
Jo 


( 49 ) 
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and proceeding exactly as before using (49) in place of (29), we immediately arrive 
at the following remarkable estimate: 


Corollary 1. For k > 1, and either 1 < a < 2 and d = 2 or 0 < a < 2 and d > 3, 
the eigenvalues {Ay of (47) defined on Q C satisfy 


k 


Y X i ^ ( 4?r ) S 

j =i 


cr d 
a + d 


r(i + l) 

\n\ 


fc 1+ t 


+ 


a a |Q |2 d r ( 

i + f) 

a-1 

d 

2(ct + d) (47r) 2 - ^ 

5 aa 1 — T ' 

J(fi)5 

(1 + C 

> d 

16 (a + d) (47T) 1 — 

aa |fi| 5 r 

I J(fi) 

(1+1 

) d /.l+s^ 

16(ct + d) ( 47 r)® — 

il(fi) 

3 lyj 

2 


(50) 


Note that this estimate also improves the main result of [19], which is simply 
the multiple of the lower bound stated in (9) [33] by cr. 
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